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ABSTRACT 

Neutrino mass matrix generated by the Zee (radiative) mechanism has zero (in general, 
small) diagonal elements and a natural hierarchy of the nondiagonal elements. It can be 
considered as an alternative (with strong predictive power) to the matrices generated by 
the see-saw mechanism. The propagation in medium of the neutrinos with the Zee-mass 
matrix is studied. The flavor neutrino transitions are described analytically. In the phys- 
ically interesting cases the probabilities of transitions as functions of neutrino energy can 
be represented as two-neutrino probabilities modulated by the effect of vacuum oscillations 
related to the small mass splitting. Possible applications of the results to the solar, super- 
nova, atmospheric and relic neutrinos are discussed. A set of the predictions is found which 
could allow to identify the Zee-mass matrix and therefore the corresponding mechanism of 
mass generation. 
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1 Introduction 



There are a number of mechanisms of the neutrino mass generation which relate naturally the 
smallness of the masses to the neutrality of neutrinos. Those include the see-saw mechanism 
0,0 different radiative mechanisms 0, (3|, ||, ||, 0, ||, the tree level mass generation by 
Higgs triplet having small induced VEV ||, mass generation by nonrenormalizable Planck 



scale interaction [10]. In general several mechanisms contribute to the neutrino masses 
simultaneously, moreover their contributions can be comparable or one of them dominates. 
Different mechanisms imply different symmetries and particle contents of the theory. 

There are some indications of the existence of tiny neutrino masses and lepton mixing 
related to the solar, atmospheric and relic neutrinos. The suggested mechanisms can rather 
easily generate the neutrino masses (mass squared differences) and mixings in some or 
even all regions of these "positive indications". What can be learned about the origin of 
neutrino masses and mixings from the existing, and mainly, from future experiments with 
solar, atmospheric, supernova neutrinos as well as from the accelerator experiments? Can 
one conclude something on the mechanism of the neutrino mass generation? In principle, 
different mechanisms result in different structures of mass matrices and consequently in 
different features of the neutrino propagation in vacuum and matter. For the identification 
of the mechanism it is crucial to study the effects of mixing of all three neutrinos. The 
propagation of three neutrinos with strong hierarchy of masses and mixings typical for 
simplest versions of the see-saw mechanism has been widely discussed before O] . 
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In this paper we will study the properties of propagation in medium of the neutrinos 
uj = (z/ e , z/ M , v T ) having the Major ana mass matrix: 



(1) 



Its crucial features are: zero (in general small) diagonal elements and strong hierarchy of 
the nondiagonal elements. Such a matrix appears in a class of models with radiative mass 
generation. The first and the simplest version which naturally results in structure (1) had 
been suggested long time ago by Zee || and we will call matrix (1) Zee-mass matrix. 

The Zee mechanism implies the existence of charged scalar field, S, singlet of the SU(2) 
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and two doublets of Higgs bosons $ M , $> d - The interactions responsible for the neutrino 
mass generation are 

Tfl — 

J2 f a ^l L Cr^ pL S^ + £ — + M ud $ u T 2 $ d S^ + h.c, 

a/3 a V d 

(a,P = e, [A, r), where ^/ql is the lepton doublet, l a R is the right handed component of the 
charge lepton, m a is its mass, and M ud is the mass parameter. As a consequence of the gauge 
symmetry the couplings of the singlet 5* to lepton doublets are antisymmetric: f a p = —fp a - 
Only one doublet, Q d , with vacuum expectation v d gives masses to the charged leptons. The 



indicated interactions generate at one loop level the elements of the mass matrix (1) p| |]12 



m ^AL( m 2 _ m 2 ) I _ ^ M ud v Ml 



where M 1; M 2 are the masses of charged scalars, v = \jv\ + v d is the electroweak scale and 
f M is the vacuum expectation of From (2) one finds the relations between the masses: 

i.e. the element m eA1 turns out to be naturally suppressed in comparison with the two other 
elements thus reproducing the structure (1), unless f a /3 have strong inverse flavor hierarchy. 

The matrix (1) appears also in some modifications of the Zee- mechanism. Instead of 
one scalar, three charged scalars, S 1 ^, S eT , S^, are introduced in that carry double lepton 
charges. Moreover, the lepton number can be violated spontaneously by the VEV of new 
neutral singlet which has the coupling Q u Q d S + S° ||. Three boson coupling can be generated 
in one loop , being suppressed at tree level by a discrete symmetry. The general properties 
of models which result in the matrix (1) are (i) the existence of one (or several) charged 
scalar fields singlets of SU2, (ii) the generation of masses of the charged leptons at tree level 
by only one Higgs doublet. Last feature explains simultaneously a suppression of flavor 
changing transitions and can be related to a certain discrete symmetry. 

As can be shown, the properties of propagation are practically the same for a more 
general form of mass matrix with nonzero diagonal elements provided such elements do not 
exceed appreciably m efl . The appearance of these diagonal elements can be related, e.g. to 
the violation of the discrete symmetry in the interactions of Higgs scalars. (The restrictions 
on flavor changing neutral currents admit rn dia9 ~ Tn efl ). They can result also from the 
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see-saw contributions or from the radiative effects related to the R-parity breaking. Note 
that such a matrix as well as the original Zee-matrix have an approximate L e + — L T 
lepton number conservation. This allows for a large mixing of e- and \i- flavors. The r flavor 
turns out to be singled out. 

The properties of the Zee mass matrix in vacuum have been studied by Wolfenstein [|13| . 
Here we consider the matter effects and confront the results with the existing data. The 
paper is organized as follows. In sect. 2 we summarize the properties of the Zee matrix and 
the oscillations induced by this matrix in vacuum. In sect. 3 the properties of the Zee matrix 
in matter will be considered. Sect. 4 is devoted to the conversion induced by the Zee matrix 
in matter. The application will be discussed in sect. 5. Sect. 6 contains the outlook and the 
discussion. 



2 Zee mass matrix in vacuum 



Instead of m^- it is convenient to introduce basic mass scale, m , e/i- mixing angle, 9, and 
small parameter, e, asQ 



m = \jm\ T + m^ r , sin 9 



fe 



m 



m 



m. 



m. 



P + P ' 

J er 1 J [it 



e/Lt 



f er f ut 



' er 1 J [it 

In terms of these parameters the mass matrix (1) can be rewritten as 



(4) 



M = m 



e sin 9 
e cos 9 
sin 9 cos 9 



(5) 



Let us summarize the properties of this matrix. The eigenvalues of (5) equal to JL3 



mi = — moesin26 l , 7713,2 = mo(±l — -esin2#), 



(6) 



1 In contrast with Wolfenstein |13], we use sm8 «-> cos 6, e — > S. 
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and the mixing matrix in vacuum which relates the flavor and the mass eigenstates Vf = Squ 
{u = U!, u 2 , U3) is 

( ^2c s s \ 

— \/2s c c 
{ l-l) 
here (and everywhere below) s = sin 9, c = cos 9. 



+ 0(e), (7) 



Since e ■€ 1, two states, v 2 and z/ 3 , are approximately degenerate, and their masses 
(~ mo) are much larger than the mass of v\. Moreover v 2 and ^3 have opposite CP-parities 
thus forming the pseudo Dirac neutrino. Their mass squared difference is 

Aml 3 = 2e sin 29m\ -C m 2 ,. (8) 

Mixing of the electron and the muon neutrinos is determined by 9, it becomes maximal at 
m eT ~ m Mr (f eT = fur)- The tau neutrino is mixed almost maximally with the combination 
v\ = sin9i' e + cos 6^; the deviation from maximal mixing is proportional to e 2 : 

1 - sin 2 2# r2 = ie 2 sin2#. (9) 

The mass squared differences for v\ component equal Am 2 2 ~ Am 2 3 ~ wig and the ratio 
of mass differences (see (6)) is determined by e: 

= do) 

Consequently, the Zee-matrix has two crucial features: it gives naturally two different scales 
for the mass squared differences and practically maximal mixing between the two heaviest 
components. 

In vacuum the propagation of neutrinos having the Zee-mass matrix results in superpo- 
sition of oscillations with two different oscillation lengths; the ratio of lengths is determined 
by (10) |13|]. For small distances, L 47r£ , /Am 2 3 , the task is reduced to two neutrino 
oscillations u e «-> z/ M with depth sin 2 29. The probabilities of 3z/-oscillations can be immedi- 
ately found from (6) (7). For example, the averaged survival probabilities P[y e — > v e ) and 
P{v„ -> Vft ) are @ 

P(u e ^ u e ) = c 4 + l -s\ P(v, - vy) = s 4 + ic 4 . (11) 

Note that the effect of the third neutrino is reduced to the factor of 1/2 in (11) which comes 
from the lost of the coherence in the maximally mixed neutrino state. 
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3 Properties of Zee mass matrix in matter: levels and 
mixing 



Evolution equation. Effective Hamiltonian. In matter the evolution equation of the neutrinos 



[y e , u^, v T ) can be written as |14 



■ dU f IT 



H — cyj? ^ matter t 



(12) 



where H is the effective Hamiltonian for ultrarelativistic neutrino, M 2 is the vacuum mass 
matrix (5) squared, E is the energy of neutrino and H matter is the matrix describing matter 
effect: H mat ter = diag[V2G Fn e , 0, 0]. Here Gf is the Fermi constant and n e is the 
concentration of electrons. In (12) we have neglected the high order electroweak effects which 



lead also to the splitting of — v T levels [17| and suggest that the concentration of neutrinos 
in medium is small (so that neutrino - neutrino scattering which in particular, generates the 
nondiagonal elements of H m attefM^ can be neglected). Substituting the explicit expressions 



for M 2 and H matter in (12) one gets 



2E 



s 2 + e 2 + p 



sc 
ec 



sc 



ec 



c 2 + e 2 es 



es 



1 



P 



\[2G pn e — =■ . 



(13) 



Diagonal elements of H determine the flavor energy levels, and as the density changes, 
there are two crossings of flavor levels. The resonance (crossing) densities are: 



p efl = cos 29, 



cos 2 9, 



(14) 



for v e — and u e — v T levels correspondently (see fig. 1). Evidently p er > p efJj ; the distance 
between crossing points (p er — p eA1 = s 2 ) is always smaller than the width of the e — /i 
resonance (Ap ~ tan2#). Consequently, there is a strong influence of the resonance related 
to the largest e — p- mixing on that stipulated by e — r-mixing (fig. 1). 



Mixing matrix in matter, level crossing scheme. Let us introduce the neutrino eigenstates 
in matter, u im , and corresponding eigenvalues Hi, (i = 1, 2, 3) as the eigenstates and the 
eigenvalues of H . The eigenstates and the flavor states are related by mixing matrix in 
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matter, S(p): Uf = S{p)v m . The matrix S(p) and the eigenvalues Hi are determined by the 
diagonalization condition S(pY H{p)S{p) = diag{Hi, H 2 , H 3 }. 

Let us find S(p) and H(p). As follows from (13) in zero order over e {m e(l = 0) the state 
v T decouples and the task is reduced to two neutrino case. The Hamiltonian H° = H(e = 0) 
is diagonalized by the rotation 



cos 9 m sin 9 m 
— sin 9 m cos 9 m 
1 



(15) 



where 9 m is the ep - mixing angle in matter fixed by: 

sin 29 

tan 20 Tl 



cos 29 — p 



The eigenvalues of H° which correspond to the rotated states = (vi m , v\ m , v T 



i • m o 
written m the units — — as 
2E 



1 + p ± ^(1 + pf - Ac 2 p 



(16) 
can be 

(17) 



The level splitting, Aif 12 = H% — H® = \J (cos 29 — p) 2 + sin 2 29, is minimal in resonance 
(14): A# 12 (p eM ) = sin2#. 

At p = 0, one gets from (17): H® = 0, and = H T = 1, i.e. the levels and H T cross 
at zero density (fig.l). Due to the strong influence of ep- mixing on er-mixing the resonance 
for v T is shifted from p er to p = 0. 

The elements of the Hamiltonian (13) proportional to e give the corrections to the above 
level scheme which become important when p — > 0, in particular they induce the mixing of 
v T — v\ m states. Using (13) and the mixing matrix (15) one finds the Hamiltonian in the 
basis of the states = (fi^i^)' Performing then an additional rotation, Si 3 (a), in 
v\ m , v T space, on the angle a determined by 



tan 2a = 



2e cos(6» + 9 m ) 
l-H°-e 2 ' 



(18) 



[u^ — > i/ ra ,^ = Si3(a)u' m , and v' m = [y\ m ,v\ m ,v' T )) one gets the Hamiltonian which de- 
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scribes the conversion at small densities 

Hi -8 

H~ # 2 ° esin(0 + m )cosa 

esin(0 + m ) cos a 1 + e 2 + <5 

Here 



+ 0(e 3 ). 



(19) 



<J = -(1 - sin 2 a + e cos(0 + m ) sin 2a 
is of the order of e 2 , the values of 9 m , H®, if° are determined in (16 - 17). 

According to (18) the angle a being also of the order of e increases with density: 



tan 2a 



p = 

p = cos 20 

p > 1 



2e cos 20, 
2e(c + s)-\ 
{ 2e/s, 

There is no resonant enhancement of a. 

The Hamiltonian (19) is diagonalized by the rotation S 2 3((3) in the {y\ m , z/)-space by 
the angle /3 determined from 

2esin(0 + 9 m ) cos a 2esin20 4e 



tan 2/? 



l-HS + 6 



(20) 



,' 2 - i 1 — e 2 sin 2 20 — ps 2 p tan 

Here we have taken into account that in the denominator 5 can be neglected everywhere 
except small region around p and that ~ 1 +ps 2 for e 2 <C p <C 1 as follows from (17). 
Obviously, tan 2(3 <C 1 for p 3> 4e/ tan0. If p 3> 1 one has ~ p and the angle /3 is even 
more strongly suppressed tan 2/5 rs 2e/p. At p = second expression in (20) reproduces the 
result (9). 

The total mixing matrix in matter that diagonalizes the original Hamiltonian (13) is 
S = 5*i2 (#m) • Sis(a) ■ S23 {(3). The eigenstates H 2 and H 3 can be found as the eigenstates 
of 2 x 2 submatrix of (19), H 1 = — 5. Two angles, 9 m and (3, undergo the resonant 
enhancement in different density regions. 

Two density regions. It is possible to divide the whole density region into two parts so that 
in each part the three level mixing is reduced to two level mixing. Indeed, let us define the 
density 

Pb = — (21) 



tan0 
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which fixes the width of the v T 



'2m 



resonant layer (Ap# = 2pb) (for p < pb one has 
sin 2 2/5 < 1/2). If pb <C 1, the regions of small and large densities can be introduced. 

(i) . Large density region: p ^> p&. Here a,j3 ~ e and in the lowest approximation 
the mixing matrix is S ~ Si2(#m)- The state i/ T decouples; the dynamics of propagation is 
determined by change of 9 m . 

(ii) . Small density region: p < p&. The change of mixing is determined by the angle /3, 
whereas two other angles vary weakly (even if 9 = 45°) coinciding practically with vacuum 

6, tan 2a ~ 2ecos26 l . Consequently, in the first approximation one has: 



values: 9 m 

S = Su (9) ■ (/£?); v \ m decouples and the dynamics of level crossing is determined by 2 x 2 
submatrix of (19). 

For extremely large densities, p>l, : 9 m ~ 7i~/2, (3 ~ and a ~ e/s so that 



S m (p » 1) 



1 
-1 ^ 

s 

e - 1 



(22) 



i.e. v e state decouples being composed of the eigenstate z/ 2m , whereas and z/ T are mixed 
with the angle a ~ e/s. 



m^, and therefore, s = c 



Other possibilities. In special case of equal / Q/ g, one has m e 

— =. The z/ e — level crossing takes place at p — 0. As before — z/2m resonance is also 
V2 

at p ~ 0. The Hamiltonian of v T — z/ 2m system is simplified: v\ m decouples since at small 
densities a ~ 0. 

Another possibility, m er > m^, is realized in case of inverse hierarchy of couplings, 

7T 

fer > ffir- Now 9 > — and v e — level crossing is at negative p; e/i-resonance conversion 
takes place in the antineutrino channel. Crossing of the v T - z/ 2m - levels occurs at p ~ as 
before. However these possibilities are disfavored by data from SN1987A WE . 



4 Neutrino transitions in medium with variable den- 
sity. 

Adiabatic conversion. Mixing matrix and the eigenvalues obtained in sect. 3 allow to get 



S 



the probabilities for oscillations in uniform medium as well as for the adiabatic conversion 
in medium with varying density. The probability of the adiabatic conversion P{y% — > Vj) 
averaged over oscillations is determined by the mixing in the initial moment, S(po), and in 
the final moment, S(pf): 

P a %-^) = £lMP/)^(Po)| 2 . (23) 
k 

Suppose v e is produced at the density po and propagates adiabatically to p = 0, then 
substituting S(po) = S^iO^) and S(pf) = So in (23) one finds the v e — > u e survival proba- 
bility 

P ad (u e -> u e ) = (1 - I* 2 ) cos 2 61 + \s\ (24) 
where 6^ = 9 m (Po) and the mixing angle 9 m is determined in (16). For p ^> 1 the mixing 

7T 1 

angle is 9^ w — , and from (24) one gets P |s 2 ; it differs by factor - from 2z/ conversion 

probability due to maximal mixing oscillation (conversion) between two nearly degenerated 

states. When p — > (0 m — > 0), the probability (24) reduces to averaged vacuum probability 

2 

(11). Note that at s 2 = - the probability does not depend on matter effects: P ad = 1/3. 
2 / 2\ 

At s 2 < - f s 2 > —J the probability P ad {y e — > z/ e ) decreases (increases) with p increase. 

The (z/ e — > z/ M )-transition probability, P ad (v e — > z/ M ), can be found from (24) by the 
interchange s 2 <-> c 2 ; the interchange c m <-> s m gives the probability of the transition 
P ad {y^ — > z/ e ). The antineutrino transitions are also described by (24); the mixing angle of 
antineutrinos in matter is smaller than that in vacuum, e.g. for p ^ 1 one has 9 m ~ and 
P(p e — > P e ) PS cos 2 0. 

Adiabaticity violation. Let us consider the general case taking into account the effects of 
adiabaticity violation and oscillations. The task is essentially simplified due to the existence 
of two different scales of Am 2 . As it was mentioned in sect. 3 in a given density region only 
one mixing angle changes appreciably, whereas two others are "frozen". Consequently, the 
three neutrino task is reduced to two neutrino tasks. In this case one can introduce partial 
adiabaticity parameter Kij that determine the probability of a jump between two given levels 
Hi, Hj as 

.{p) l3 = AH - {P) ■ < , (25) 
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where ip = 4^ determines the change of the level mixing in a given density region and 
AHij(p) is the level splitting. 

In the region of large densities, p ^> 4e/ tan#, the change of mixing is stipulated mainly 
by 9 m , i.e. ip = 9 m , and in the e — p resonant point one gets using (16): 

_ 2sin 2 2fl mf 

As can be shown for a not too small mixing angle 9 the adiabaticity for the 2-3 levels is 
fulfilled much better than for the 1-2 levels so that with increasing p the adiabaticity starts 
to be broken first for the 1-2 levels and then for the 2-3 levels. In the case of complete 
adiabaticity for all levels the neutrino state produced at p 3> 1 as v e follows the H 2 level 
(fig.l). If the adiabaticity of v% m — v T system is broken then neutrino evolves along the TT\ 
trajectory which is actually very close to H 2 for p ^> e/s 2 . Consequently, in this region it 
does not matter whether the 2-3 level adiabaticity is broken or not. In case of a strong 
adiabaticity violation for 1-2 levels in ep-crossing region the neutrino state follows the 
z/ e -trajectory. 

In the region of small densities, p < 4e/ tan#, one has ip — f3. Here 9 m 9, a « 0. The 
adiabaticity parameter for v\ m — v T levels in resonance (p ~ 0): 

R 32e 2 cos 2 9 ml 

is much smaller than K\ 2 {p = 0), i.e. the adiabaticity can be broken for v\ m — v T levels, 
whereas v lm propagates adiabatically. 

Another circumstance which simplifies the task is the maximal mixing of v\ m and v T 
levels at p = 0. This ensures that the probabilities of the transitions with zero final density, 
Pf = 0, averaged over the oscillations do not depend on the adiabaticity condition in the 
u 2m ~ u t system (see Appendix). 

The probabilities of conversion. Keeping in mind possible applications to the solar, su- 
pernova and the atmospheric neutrinos we will consider the propagation of the electron 
neutrino, u e , produced at some density po towards zero density, pf = 0. In general the adi- 
abaticity may be broken in ep- resonance region as well as in the region of small densities 
(p ~ 0). Leaving the medium the neutrinos will oscillate in vacuum, and moreover, the 
oscillations induced by small mass splitting may not be averaged out. 
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According to (15) the decomposition of the initial neutrino state over the instantaneous 
eigenstates is 

i/ = v e ss cos 6° m is° lm + sin 6° m v° 2m , (28) 

where 6^ is the mixing angle in the production point (the admixture of u 3m , being of the 
order of e, is practically unessential). Let us introduce some density p' (p& < p' <C 1) so that 
at p < p' the 1-2 level adiabaticity is restored or the change of the 1-2 mixing is negligibly 
small. As the result of propagation over the large density region one gets then at p' the 
state 

u(p') = (A u c° m + A 21 s°Je^uUp) + (A2C° m + A 22 s°M m (p'), (29) 

where Ay (i,j = 1, 2) are the amplitudes of transitions between the levels v\ m and v\ m , 0'm 
is the phase, (9 m = cos 9^ etc.. The amplitudes satisfy the relation: 

\A 12 \ 2 = \A 21 \ 2 = 1 - \A 22 \ 2 = P 12 . (30) 

The jump probability P\ 2 can be approximated by the Landau-Zener probability [0| (or its 
modifications): 

Pi2 ~ Plz = ezp(-~K&)> (31) 

where k^ 2 is the adiabaticity parameter in resonance (26). For the adiabatic transitions one 
has Aij = 3{j • 

In the region of small densities, p < p', the state v lm f=3 v\ m propagates adiabatically so 
that its admixture does not change, and at zero density one gets v\ m P3 v\. The evolution 
of v 2m state is described by matrix (19). The result of its propagation to zero density can 
be presented as 

4m ^ |a 22 |e^ 2 + |a 23 h , (32) 

where a 22 and a 2 3 are the amplitudes of the transitions v 2m v 2 and v 2m — > u 3 correspon- 
dently. In case of the adiabatic propagation a 22 = 1, a 23 = 0. Note that neutrino crosses 
only half of the 2-3 resonance region and therefore the jump probability equals half of Plz'- 
P 23 = | a-23 1 2 ~ \Plz- Combining (32) and (29) one gets the neutrino state at the edge of 
medium: 

i/(0) = (A llC ° m + A 21 s° m ) e l ^v x + (A 12 c° m + A 22 s° m ) (\a 22 \e^u 2 + \a 23 \u 3 ) . (33) 



11 



Further propagation in vacuum results in changes of phases onlyQ: <f> m — > = <j> m + mc and 

Am 



' ' - C + &,ar, where 



III ' III r/ir- ^ 

23 L (34) 



2E 

is phase difference acquired at a distance L in vacuum. Using (33) and the vacuum mixing 
matrix (7) one can obtain the probabilities of different transitions. In particular, [v e — > v e )- 
survival probability averaged over short length oscillations (phase 0') is 

P[v e ^v e ) = \< u e \u(0) > | 2 = c 2 (cos 2 e° m - P 12 cos 29°J + 
+ |s 2 (sin 2 ^ + P 12 cos2# m ) (1 + 2\a 22 a 23 \ cos0) , 

where we have taken into account (30) as well as similar normalization condition for Oy. 
The probability averaged over 0: 



P(l>. — K ) fl-^ 2 



cos 2 ^ -P 12 cos2# ! 



-,2 



+ \ (36) 



does not depend on ay at all in accordance with the general statement proved in the Ap- 
pendix. At P\ 2 = (adiabatic propagation in the region of large densities) the result (36) 
coincides with (24). Using (36) we can rewrite the probability (35) as 



P[v e ^u e ) = P + s 2 (sin 2 6° m + P l2 cos 26° m ) ^ P 2Z (l - P 23 ) cos 0, (37) 

where P 23 = 1 0-23 1 2 - If the adiabaticity for v\ m — v T levels is strongly broken then one has 
«22 = a 23 = 75 or -^23 = \ which corresponds to z/° m — > = u ' 2 ^ 3 transition in medium. 
Substituting P 23 = | in (37) one gets 

P{y e -> v e ) = P + ^s 2 (sin 2 ^ + P 12 cos 2^) cos 0. (38) 



2 

In this case the depth of vacuum oscillations is maximal. 

Energy dependence of the suppression factors. Let us consider the dependence of the prob- 
ability (38) on the neutrino energy (fig.2). P(P) is the oscillating curve inscribed into 
the band between p max and P mm (the oscillations are stipulated by change of 0). For 
cos(4r|r-L) = 1 one gets from (38) 

P max = -+(-- P 12 ) cos 29° m cos 29 (39) 



2 V2 



2 The lost of coherence due to wave packet spread is reduced to the averaging effect. 



12 



which coincides with 2i/-probability P<iv, i.e. P 2u gives the upper bound for 3z/-survival 
probability: P 3u < P% v . The width of the band, 

A p _ pmax _ pmin = ^ + p^ ^ ^Oj^ ^ 

is proportional to sin 2 9, and consequently, with diminishing 9 the 3^-probability converges 
to 2 ^-probability. For neutrinos propagating in matter with monotonously changing density 
pmax _ _p 2i/ (_g) hag form of pit (fig.2) |TJ]]. Outside the pit the probability approaches 
1 on the right hand side and the vacuum value, sin 2 2#/2, on the left hand side. The 
position of the left (adiabatic) edge of the pit, E aa -, is fixed via the resonant condition by 
the density in the production point. The position of the right (nonadiabatic) edge E na is 
determined by the adiabaticity condition. As follows from (40), AP = s 2 sin 9° m in the 
adiabatic region. In particular, for small energies outside the pit, where the matter effect 
is weak, one has AP = s A . Maximal width of the strip, AP = s 2 , is at the bottom of the 
pit when po >> Pr, therefore at the bottom: p mm = 0. In the nonadiabatic region AP 
decreases with enhancement of the adiabaticity violation: AP ~ s 2 (l — P12). 

The position of the first (broadest) minimum of the oscillating curve, E m , is determined 
from the condition — — L, where l v is the oscillation length in vacuum. Explicitly one 
A771 2 L 

has E m = — — ; the first maximum is at E m /2 etc.. Mutual position of the pit and the 

2ir 

modulating curve is fixed by the ratio of energies 

E ad = 27rcos2fl Am 2 12 _ Am 2 2 / 
E m ~ LV2G F n Am 2 23 ~ Am| 3 L' 

where Iq = 27r/v / 2GVn cl is the refraction length in the neutrino production point. If E rn ^> 
E na > E ad , fast oscillations of the modulating curve are practically averaged in the energy 
region of the pit and P ~ P. For E m ~ E ad one predicts the observable modulations of the 
pit. If E m E a long length oscillations are not developed and P — * P 2u - 

Similarly one can analyze the transitions of and v T . In particular, v T is converted 
mainly in the region of small densities, and the probability of v T — > v e transitions averaged 
over large scale splitting equals 

.2 r 



S 



1 + a/P 2 3(1-P 2 3)cos. 



(42) 



For the antineutrino channel (negative p) there is no v e — level crossing; one can consider 
adiabatic evolution of v\ m and v^ m states and as the result P{y e — > z/ e ) « P ad = cos 2 9. 
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If f eT ~ f/j/r: the e/z-mixing becomes maximal; both resonant regions are at p = 0. The 
above consideration (reduction to two neutrino tasks) is valid due to difference in the widths 
of the resonance layers. Indeed in v T — v\ m crossing region the change of the v e — mixing 
is negligibly small. 



5 Applications 

The results of the solar neutrino experiments can be reconciled with predictions of the 
Standard Solar Models in terms of the resonant flavor conversion v e — > v^iv^) with pa- 
rameters (two neutrino mixing): Am 2 = (0.4 — 1.2) • 10~ 5 eV 2 , sin 2 26 = (0.1 — 1.5) ■ 10~ 2 
"small mixing solution" and Am 2 = (0.5 - 3) • 10" 5 eV 2 , sin 2 26 = (0.60 - 0.85) "large 



mixing solution" [20 1 . However taking into account the uncertainties of both the predic- 
tions and the experimental data, one should consider more wide region of the parame- 
ters. The deficit of in the atmospheric neutrino flux testifies for oscillations z/ M — > i/ T , 
with Am 2 = (0.5 - 3) ■ 10~ 2 eV 2 , sin 2 26 = (0.4 - 0.6) or for oscillations ^ -> u e with 
Am 2 = (0.5 - 3) ■ 10~ 2 eV 2 , sin 2 26 = (0.3 - 0.8) [gl|]. The formation of large scale structure 
of the Universe implies the existence of hot component of dark matter and the relic neutrinos 
with m ~ (2 — 7) eV can play such a role p2| . 



Let us confront these results with predictions of the Zee-model for different values of 
m . As follows from (2, 4), basic mass scale is m = 7 ■ 10 4 ■ {iy/fl, + fl T ) eV, therefore 
depending on values of couplings fij as well as the parameters of the scalar sector one can 
get for m any value below ~ 10 4 eV. For f efM < f er < f^ T the ratio of the mass squared 
difference (8) may be in the range 10~ 5 — 5-10 -3 , unless one introduces very strong hierarchy 
of couplings. Also under this condition the mixing angle is in the region sin 2 26 = 10~ 3 — 1. 

Three regions of tuq values are of special interest. 

Cosmologically interesting mass scale: mo ~ (1 — 30) eV. The components i>2 and v% with 
masses mo, can form the hot dark matter. Since Amf 3 ~ Am 2 2 ~ m\ > leV 2 , the mixing 
angle, 6, is restricted by the accelerator oscillation experiments: e.g. sin 2 26 < 2 ■ 10~ 3 at 
m = 3 eV, sin 2 26 < 6 • 10~ 3 at m = 1 eV etc. [23|. For sin 2 26 > 10~ 5 one expects strong 



resonant conversions v e — > z/„ and — > v e in the inner parts of the collapsing star. Such 
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a conversion results in permutation of the v e - and v^- energy spectra and therefore in the 
increase of average energy of the electron neutrinos. This will have two consequences: (i) the 
increase of the energy release due u e — e- scattering which may help to expel the envelope, (ii) 
the formation of the proton-rich medium due to dominant v e n — > ep scattering. The latter 
will forbid the r-processes responsible for nucleosynthesis of heavy elements EM . If the inner 



part of collapsing stars is the only place of the r-processes, then the indicated conversion 
should be suppressed and one gets the bound on mixing angle sin 2 29 < 10~ 4 — 10 -5 [E4 



The value of Am| 3 , can be naturally in the region responsible for the atmospheric neutrino 
problem. The suppression of the muon neutrino flux due to — z/ T -oscillations is determined 
by the averaged vacuum probability (11) and taking into account the indicated bounds on 
9 one gets: P ~ 0.5 which is actually outside the region of the best fit of all the data. 
The suppression of the solar z/ e -flux fixed by averaged vacuum probability (11) is very weak: 

P{v e U e ) > 0.92. 

Atmospheric neutrino mass scale: niQ ~ (0.3 — 1) • 10 _1 eV. j«q is in the region of the solution 
of the atmospheric neutrino problem Am^ pa m 2 , ~ (10~ 3 — 10 _2 )eV 2 . The deficit of can 
be explained by «-> v e - oscillations, with sin 2 29 = 0.3 — 0.8. 

For A77I23 = (10~ 6 — 10~ 4 ) eV 2 the suppression factor for solar neutrinos is determined by 
the adiabatic probability (24) which coincides at 9 m pa 9 with averaged vacuum probability 
(11). Solar neutrino spectrum is outside the pit at small energies. Although Am| 3 is in 
the region of strong matter effect, the averaged probability P , practically does not depend 
on matter density and on the neutrino energy as well. For values of 9 needed to solve the 
atmospheric neutrino problem (at s 2 < c 2 ) one gets from (11) P = 0.56 — 0.83. Taking into 
account an additional contribution to the z/e-scattering from neutral currents one predicts 
the following suppression factors (R = data/ SSM) for Ga- , Ar-, production rates and 
z/e-signal: 

R Ge = R Ar ~ P ~ 0.56 - 0.83, R ue ~ 0.62 - 0.86. (43) 

There is no distortion of energy spectrum of boron neutrinos. The predictions (43) fit rather 
well all the results except the one of Homestake experiment. 

In case of very small splitting Am 2 3 <ti 10 -9 eV 2 , there is no averaging over long 
length oscillations and the probability is modulated with the amplitude s 4 /2 (see (40) and 
further discussion). For the indicated mixing angles the amplitude of modulations is ~ 
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(0.25 — 4.5) • 10~ 2 . So, one may expect up to ~ 10% distortion of boron neutrino spectrum 
and change of the relations (43). For example, at sin 2 29 = 0.8 and for certain values of 
Am 2 one may get Rc e = 0.54, Ra t = 0.51, and R Ue = 0.58. Further increase of sin 2 29, 
although enhances the amplitude of modulations, results in the stronger suppression of the 
atmospheric z/ M -signal as well as the signal in gallium experiments. 

As in the previous case one predicts strong resonant conversions v e — > and — > u e 
of the neutrinos from the collapsing cores of stars. The corresponding permutation factor 
which characterizes the interchange of v e - and v^- energy spectra [T5[ equals p = 0.75 — 0.90. 
Moreover, since the mixing is rather large one expects an appreciable permutation of the 
antineutrino spectra: v e <-> v^; the permutation factor p = 0.08 — 0.25 results in even better 



description of data from SN1987A [27|, [15 



Note that for indicated values of mo the conversions take place now in more external 
layers, so that there is no problem with r-processes. 

Another value of the mixing angle which follows from explanation of the atmospheric 
neutrino data (s 2 > c 2 , this corresponds to the consideration in |25fl ) gives P = 0.17 — 



0.44 and P ue = 0.29 — 0.51. The result contradicts to the observed gallium production 
rate. Moreover the e — /i-resonance is in the antineutrino channel which results in strong 
permutation ( p = 0.6 — 0.8) of the u e — , v^— energy spectra. The SN1987A data give the 
bound p < 0.4 



Solar neutrinos mass scale: m — (10~ 2 ~ 10~ 3 ) eV. the largest mass splitting, Am 2 3 ~ 
(10~ 6 — 10~ 4 ) eV 2 , is in the region of the resonant effect inside the Sun. If Am 2 3 >> 
10 _10 eV 2 , averaging over long length vacuum oscillations takes place and the suppression 
factor for z/ e -ffux is determined by the probability P (36). It can be rewritten as: 

r sm 2 9 m + P 12 cos29 m , (44) 



where P 2v is the suppression factor for two neutrino mixing. One can easily construct the 
3^-suppression pit using (44) and the results for two neutrino mixing. For small 9 the effect 
of the third neutrino is negligibly small and P ~ P 2u . The Zee-mass matrix reproduces 
small mixing solution of the solar neutrino problem for two neutrinos. The deviation from 
the 2z/ case is of the order of 10~ 3 and to distinguish the Zee mechanism from the other 
mechanisms one can take into account the following facts. In the considered case there 
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is no oscillation solution of the atmospheric neutrino problem as well as no appreciable 
contribution of neutrinos to hot dark matter. Moreover, there is no manifestation of the 
third neutrino in the experiments with supernova neutrinos. 

Large mixing solution is absent unless one admits large original flux of boron neutrinos. 
Indeed, since at the bottom of the pit P = sin 2 9/2 one needs two times bigger value of sin 2 9 
in comparison with 2v case to get the same suppression of Kamiokande signal, e.g. instead 
of sin 2 29 = 0.7 one should take sin 2 29 = 0.994. However in this case the fluxes of the low 
energy neutrinos are strongly suppressed: in the pp-neutrino region one gets P ~ 0.375. 

The situation is essentially different if Am 2 3 < 10~ 9 eV 2 , and consequently there is 
no averaging over the long-length oscillation at least in some part of the suppression pit 
(fig.2). One of the most interesting configurations is shown in fig.3 which corresponds to 
Am^/Am 2 ^ ~ 10~ 5 and sin 2 29 = 0.8 — 0.9. First minimum of oscillating curve is at 
adiabatic edge, first maximum is outside the pit. If the Be-neutrino line is in first minimum 
of oscillating curve then the boron neutrinos are at the bottom of suppression pit outside 
the minimum, and the detected part of the pp- neutrino spectrum is in the first maximum. 
The signature of such a scenario is the strong suppression of the Be-neutrino flux, and the 
absence of the distortion of the high energy part of the boron neutrino spectrum in contrast 
with small mixing solution for two neutrinos. For pp-neutrinos one can get the suppression 
0.55 - 0.60, so that total Ge-production rate could be about 50 - 70 SNU in agreement with 
present data. 

Another possibility corresponds to the modulating curve shifted to larger energies, so 
that boron neutrinos are in the first maximum (pp-neutrinos are in averaging region) and 
the beryllium line is in the fastly oscillating region. If the oscillating curve is shifted to 
lower energies then one expect the distortion of the pp-neutrino spectrum. 

Obviously in this case there is no solution of the atmospheric neutrino problem. For 
the indicated values of parameters one may expect complete or partial conversion of the 
supernova neutrinos v e <-> depending on density profile of star. 

The loop diagrams, similar to those generating the neutrino masses, will generate also 
the transition magnetic moments of the neutrinos. These moments are however restricted 



by fj, < 3 ■ 10 /iB(mo/leV), where [ib is Bohr magneton f26| , so that even for mo ~ 10 eV 



the effects of spin-flip on the solar and atmospheric neutrinos are negligibly small. 
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Implications to the parameters of the Zee model. Let us comment on possible implications 
of the above results to the original Zee-model. The expression for m can be rewritten as 



Moreover, the dimensionless parameter / can be of the order 1 when all mass parameters 
entering I are at the electroweak scale and v u ~ v d . According to (45) a value 1 = 1 gives 
the lower bound on f^ T . In principle, f a p can be as large as 1; f^ T ~ 1 gives the lower bound 
on /, and consequently, the upper bound on the mass of charged scalar M 2 . Let us estimate 
the range of the parameters. 

For mo in cosmologically interesting domain one has f^ T > (0.13 — 4) • 1CT 4 and M 2 < 
(2 — 8) • 10 4 GeV. The indicated values of mass squared differences and accelerator bounds 
on 9 correspond to f efl ~ f eT ~ (0.03 — 0.1)/ MT . If mo is in the region of the atmospheric 
neutrino problem one gets f^ T > (0.4 - 1.3) • 10~ 6 and M 2 < (0.8 - 2.4) • 10 5 GeV. All 
constants can be of the same order: f^ T ~ f eT ~ f efl . Alternatively, if there is no averaging 
over long length oscillations f efl = lO -4 /^. For mo in the region of the solutions of the 
solar neutrino problem the parameters are f^ T > (0.13 — 1.3) ■ 10~ 8 and M 2 < (3 — 4) • 10 6 
GeV. Small mixing solution implies f er < 0.1 f^ T . Large mixing solution with modulations 
by long length oscillations is realized when f eT ~ / MT and / 6(U ~ 10~ 2 / MT . 

6 Discussion and Conclusions 

1. Zee- matrix in matter (for e 1) is an example of "solvable" 3z/-task. This allows to 
trace some interesting features of the dynamic of propagation, in particular, the effect of 
strong influence of one resonance on another. Dominant v e — mixing shifts the resonance 
for v T to zero density, thus changing a naive picture of level crossings. 

2. In the supersymmetric generalization of the model new diagrams appear with sleptons 
and higgsino in the loops. General structure of mass matrix is the same as (1). Note that 
Zee singlet can be embedded in the SU(5) GUT scheme by introducing the antisymmetric 
10-plet of scalar s. 

3. Zee-matrix can be considered as an alternative to the one generated by the see-saw 
mechanism. Let us note for a sake of completeness that in principle the matrix (1) can be 




(45) 
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reproduced by the see-saw mechanism too. For this one needs a special structure of the 
Majorana mass matrix of the right components, Mr. Namely, in the Dirac neutrino basis 
(where the Dirac matrix is diagonal) Mr should have zero determinants of three submatrices: 
MaMjj — Mfj = 0, (ij = 1,2,3). 

4. The mass matrices generated at the one loop level in the model with explicit R-parity 
violation |j] also differ from (1). The generic property of these matrices is the existence of 
nonzero diagonal elements which are not suppressed in comparison with nondiagonal ele- 
ments. If purely lepton couplings dominate over couplings of quark and lepton supermul- 
tiplets, the elements m er , m^ T , are naturally suppressed with respect to others by factor 
of m^/m T . Moreover, the element m TT has even stronger suppression: (m M /m r ) 2 . Such a 
matrix allows to explain both solar and the atmospheric neutrino deficits according to 



the scenario suggested in [^5 

Similar structure of mass matrix appears in the model with two loops generation of the 
neutrino masses 0. 

5. Practically all the extensions of the standard model imply the existence of neu- 
tral fermions which can play the role of right-handed neutrinos. In this case the see-saw 
mechanism is obtained and in addition to the radiative mass terms one gets the see-saw con- 



,2 



tributions to the neutrino mass matrix: m ss ~ mpM' 1 ™^. The biggest term is m TT « —ff-, 
where the Dirac mass, 17130, can be as large as the top quark mass. Depending on values of 
parameters which may be in the range rri£> = (1 — 10 2 ) GeV and M = (10 10 — 10 18 ) GeV one 
can get negligibly small, or comparable with erriQ, or even dominant see-saw contribution. 
For example, at m 3 £, ~ 100 GeV and M = 10 18 GeV one has m TT ~ 10 -5 eV which is of 
the order of em® for m < 10 -3 eV. If M < 10 16 GeV this contribution becomes dominant. 

Let us comment on the simplest possibility when only one additional element, m TT , is 
important. Such a situation is realized if the Majorana mass matrix of the RH-components 
has no strong hierarchy. In this case one can get easily both large \l — t- mixing angle, 9^, 
needed to solve the atmospheric neutrino problem and small mixing solution of the solar 
neutrino problem. Indeed, now /xr-mixing is of the order of mo/m ss and large value of # MT 
follows from the fact that the see-saw and the radiative contributions are of the same order. 
Moreover, for small efj, mixing the following relation between the ratio of masses and the 
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fir- mixing exists: 




1 - cos 26 



1 + cos 26 



(46) 



For sin 2 26^ T > 0.4 one gets ^§ > 0.01 which is roughly consistent with desirable value. The 

m 3 

ratio can be further corrected if one takes into account the see-saw contributions to other 
elements of matrix. 

In conclusion, we have considered the properties of propagation of the neutrinos with Zee- 
like mass matrix. Crucial features of the matrix are zero (small) diagonal elements and 
natural hierarchy of nondiagonal elements. It can be considered as an alternative to the 
matrix generated by the see-saw mechanism as well to other radiative mechanisms. The 
probabilities of the conversions induced by the Zee-mass matrix as the functions of the neu- 
trino energies can be represented as two neutrino probabilities modulated by the oscillating 
curve related to the long length oscillations. 

The Zee-mass matrix does not allow to explain all three problems related to the so- 
lar atmospheric and relic neutrinos simultaneously. If two heavy components are in the 
cosmologically interesting domain then one can get for atmospheric muon neutrinos the 
suppression 1/2 with no appreciable effect for solar neutrinos. The Zee- mass matrix allows 
to fit well the atmospheric neutrino data if the masses of heavy components are in the region 
~ 0.1 eV. In this case one predicts energy independent suppression of u e flux at the level 
0.6. For smaller mo the matrix can reproduce with high precision small 2z/-mixing solution 
of the solar neutrino problem without any appreciable manifestations of the third neutrino. 
For large mixing 6 new configurations of the suppression appear. In particular, one may 
have strong suppression of the beryllium line and energy independent suppression of the 
high energy part of boron neutrino spectrum. In general, the modulations of the smooth 
energy dependence of the probabilities for 2z/-case are expected. 

These features may allow to identify the Zee-mass matrix and consequently the corre- 
sponding mechanism of neutrino mass generation. 
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Appendix 

Suppose the neutrino propagates in medium with density decreasing from pi in the initial 
point to zero. Let us show that the averaged oscillation probability does not depend on the 
adiabaticity condition for the levels having maximal mixing in vacuum. 

Indeed, the initial neutrino state can be written as 

Vi = av Xm + bv 2m + cv 3m . 

Suppose u lm - state propagates adiabatically, whereas the adiabaticity for u 2m — v 3m m ay 
be broken. The latter means that there are the transitions v 2m v 3m while the neutrino 
propagates to zero density. At zero density the neutrino state will have the form 

Uf = av x + b'v 2 + cV 3 . 

In general and \c'\ differ from |6| and |c|, however the normalization condition implies 
that 

1 7 1 2 i I 1 2 1 7 / 1 2 i | / 1 2 i I 1 2 

\b\ + \c\ — \b I + \c\ = 1 — \a\ . 

The probability to find the neutrino v a = a a ui + b a v 2 + c a u 3 in final state averaged over the 
oscillations equals 

Maximal mixing of v 2 and z/ 3 in i>f means that b a = c a , and then taking into account the 
normalization condition one gets 

p = I4«l 2 + l^l 2 (l^l 2 + M 2 ) = I4«l 2 + K\\i - M 2 ) 

which does not depend on changes of b and c. 
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Figure Caption 



Fig.l Energy levels of the neutrinos with Zee- mass matrix as the functions of the 
matter density p (full lines). Dashed lines correspond to the levels in zero approximation 
over e. Dotted lines show the flavor levels, a). Small mixing angle 9, b). for 9 = 45°. 

Fig.2 Survival probability P{y e — > u e ) as the function of the neutrino energy for 
different values of mixing angle 9 (solid lines) a), sin 2 2$ = 0.64, b).sin 2 2# = 0.84. Dashed 
lines correspond to P max and P mm ; the averaged probability is shown by dotted line. 

Fig. 3 The suppression factor for large mixing and nonaveraged vacuum oscillations. 
Also shown is the solar neutrino spectrum (hatched). 
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